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Abstract 

We look at the contribution of Quark-Mass-dependent twist-4 operators to 

the second moments of spin averaged structure functions and the Bjorken sum 

rule. Its contribution is non-negligible in the former case due to large Wilson 

coefficients. We also discuss the values of the twist- 4 spin-2 nucleon matrix 

element within present experimental constraints. 
PACS number: 13.60.Hb 



Typeset using REVT^ 

*0n leave of absence from Department of Physics, Yonsei University, Seoul 120-749, Korea 

1 



I. INTRODUCTION 



Recent precision spin averaged lepton-hadron Deep Inelastic Scattering(DIS) data at 
CERN and SLAC p provide us with useful information about higher twist effects 

in spin averaged structure functions(F2 and Fl). From these data, we have previously 
extracted the higher twist effect in the second moments of the structure functions and 
obtained a model independent constraints on the nucleon matrix elements of the twist-4 spin- 
2 operators. These matrix elements provide useful information about correlations between 
partons inside the nucleon. 

In addition, polarized DIS provides crucial information on the spin structure of the 
nucleon The lowest moment of this spin structure function is related to the nucleon 

matrix element of the axial current. To actually estimate the quark contribution of the 
proton spin and to understand differences in experimental measurements at different scales, 
it is necessary to estimate the higher twist-effects |T0|-[12|. 

Higher twist operators appearing in both spin averaged moments and the Bjorken sum 
rule have been considered for some time ||T3|-[16|. In this work, we study the contribution 
of twist-4 operators, that are proportional to the quark masses. These operators have been 
neglected so far because of negligible up and down quark masses compared to the hadronic 
scale. For the Bjorken sum rule, the additional twist-4 operator is (m^757Q,g) and indeed 
negligible. However, for the second moment of spin averaged structure function, the quark- 
mass dependent operator {qDaDpmq) contributes with a large Wilson coefficient compared 
to other twist-4 operators and can not be neglected. 

In section 2, we will first discuss the Bjorken sum rule. The relevant Operator Product 
Expansion (OPE) for the time ordered vector current is derived using the method based on 
the Fock-Schwinger gauge, which has been used extensively in QCD sum rule calculations 



| r?| and DIS |[T^|T^. In section 3, we will discuss second moments for the spin averaged 
structure function and discuss the importance of quark mass dependent operators in relation 
to the recent estimates for the value of total twist-4 effects. We will also discuss the values 
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of twist-4 spin-2 matrix elements within present experimental constraints. 



II. BJORKEN SUM RULE 

OPE 

Let us consider the time ordered correlation of the electromagnetic current. 

T,M = ^ I d'xe'^^T[f;-{x)jr]) (1) 



The perturbative part is symmetric in fiu and the lowest antisymmetric part contributing 
to the Bjorken sum rule is related to the axial current. A simple way to obtain the twist-4 
spin-1 contribution, is to make use of the Fock-Schwinger gauge for the external gauge fields. 



We will follow similar steps given in ref. [14 



Here, the S^^\x,0) are quark propagators in external gauge field with i being the sum of 
the dimensions of external gauge field and quark mass. These are summarized in Appendix 
A. There are only two independent twist-4 operators contributing to this order. 

Ua = qg Gap IpQ'^q, La = q'jalbm'^Q^q, (3) 



where, Gai3= \^ai3nvG^i, and conventions for 60123 and 75 follow that of Muta |T9[. Q is the 
flavor SU(2) charge matrix. To express the second and third term in eq.(0) in terms of these 
independent operators, we have to extract the spin-1 part of the following operators. 

q Da [7/3, Islmq, q DaDp l^lsq (4) 
This can be done by assuming the spin-1 part of the operators to be. 



OapS — C^aflT ^ = gafsXs + QpsYa + QasZp + iea/SSaUa- (5) 



Taking all traces and using the following identities, 



We find, 
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q Dec [lp,l5\mq = -eapsaLa 

^ ^ 5 1 1 

q DaDp iblsq = 9ap-^iL - U)s + ga5-^{U - L)p + gpsj^iU - L)„ (7) 

Using these results, the final answer including the leading order result gives, 

^ i i 16 20 

q q t) t) 

where the leading operator is, 

= Qlpl^Q'^Q- (9) 

Correction to the Bjorken sum rule 

The contribution of is easy to estimate. One should note that to a good approxima- 
tion, Lq, = m^Va, where is the average quark mass squared. So the contribution of this 
operator compared to the leading result would just be Since m is 5 ~ 10 MeV^, its 

contribution is indeed negligible down to very small Q^. 

III. SECOND MOMENTS 

OPE 

To obtain the second moment, we have to consider the OPE to dimension 6 such that 
the contribution, where one quark is connected, is given by. 
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+ ^x''x^{q DaDp i^S\x, 0)7.Q'g)] + DaDf^Ds l^.S\x, 0)7.Q'g)] 

+ (10) 
There are four independent twist-4 operators appearing to this order. 

Aap= gq[D^,Gai,hf3Q^q = g'{qiat''q){qipt''Q\) 

Cap = qDaDfsmQ'^q 

Dal3 = 9q[Da, G^f3]-ff,Q'^q (11) 

Here, the operators are assumed to be symmetric and traceless with respect to the lorentz 
indices, are the generators of SU(3) color matrix normahzed to Tr(t")^ = |. 

The last operator does not contribute to the final answer. Similar steps have to be taken 
as before to reduce eq. (0) into the above independent form. Useful identities needed in the 
intermediate stages are given in Appendix B. The final answer is, 



+e„„(i/l - ji? - ^C)l, (12) 

Here, A,B,C are the spin-independent part of the spin- averaged nucleon matrix elements 
of the twist-4 spin-2 operators in eq.(|lT]) such that (A^^)Ar = {PaPp — \M^ga/3)A^ with 
= A,A'^ = B, = C. The polarization tensors are defined as e^,^ = g^^, — q^q^/q^ and 
d^^u = -PfzPuq'^/ip ■ qf + {Pi^qu + Puqfj.)/p ■ q - Of^u with = is a 4-momentum of 

the nucleon with p^ = M^). There is another contribution to twist-4 part coming from 
a disconnected four quark operator. This part does not have any additional quark mass 
dependent term, such that the full twist-4 part in lowest order in is. 



^;r*-' = ^^..M4 + r;,, (13) 



where Ma is defined from 



M^iPaPp - \M^N9af3) = {{q9^1al^T^q){qg^lpl^Ty)) n (14) 



Eq.([l3l) was obtained without the quark mass dependent operators in ref.( 0,|T^]) 



Estimate of C 

C is related to the second moment of the structure function e(x) At present, there 
is no direct experimental way of measuring e{x). To estimate C, we will make the following 
approximation. 

{qD^DpQ^mq)M ~ -^P„"P,"(m„MM)^ - ^P^PtimJd)^ (15) 

Here, P"(P'^) is the average momentum carried by each u{d) quark inside the nucleon, which 
should be | of the nucleon momentum P such that the total momentum carried by three 
quarks is roughly |P. Assuming {rriuuu) ~ 2{mddd) in the proton and using S^tv ~ 45 ± 10 
MeV [^, we have in the covariant normalization 



{muUu)N = '2{mudd)]y = 0.06 GeV^. (16) 

This implies 

{qD^DpQ^mq)N ~ (P^P/s - \Mlg^p) x (-0.0008) GeVl proton (17) 
So C = -0.0008 GeV^ for the proton. 

Experimental Data on twist-4 effects in F2 

We have recently extracted the twist-4 effects in the second moment of the structure 
function F2 0. Assuming, 

P2,L(a:,Q^) = F:^2\x,Q') + i^F^^f (x,Q2), (is) 
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where -FJ^^ takes into account the target mass correction. We found that at = 5 GeV^ 

F^=^{x)dx = 0.005 ± 0.004GeV^ (proton) 

Jo 

= ^(M4 + ^A + l5-^C), (19) 

where the second hne shows its relation to the matrix elements including C. The quoted value 
was obtained from estimates of F2^^{x) in ref. which in turn was obtained by fitting 
the BCDMS and SLAG data on the hydrogen target. We also quoted a value for the neutron 
which was based on the analysis of NMG group for higher twist effect in F^/Ff' = 2F2/ F2—I 
0. This ratio was obtained from the measurement in the deuteron and proton targets 
without taking into account any nuclear effects in the deuteron structure function. However, 
recent analysis showed |]2^ , p5| that the nuclear effect is not negligible and accounted for 



almost all the deviation from the the perturbative ratio. So we will not quote the neutron 
value and hereafter discuss the case of the proton only. 

Using the value for G From eq.(|l^ and substituting in eq.(|l9D, we see that the contribu- 
tion from G accounts for 26 % of the experimentally estimated twist-4 effects and can not 
be neglected. 

Experimental Data on twist-4 effects in Fl 

In ref. 0, we made similar estimates for the twist-4 effect in the second moment of Fl. 
This was based on the SLAG data analysed in ref. |^. At = 5 GeV^ we found 



r Fl=*dx = C 8K^F^^{x)dx = 0.035 ± 0.012GeV^ (Proton) 
Jo Jo 

2M 8 4 ^' ^ ^ 

where the second line shows its relation to matrix elements. 

Using the value for G From eq.([T7|) and substituting in eq. (^0|) , we see that in this case, 
contribution from the quark mass dependent operator accounts for only 5%. 



An estimate of the proton matrix element of M4, A and B 



7 



Let us now estimate the value of M4, A and B for the proton. We will assume C = 
—0.0008 GeV^. There are only two constraints (eg. (p!9|) and eg. (|20|) ) and 3 parameters to 



determine. However, there is one positivity constraint ( for the four guark operators. 
The proof is simple. Suppose the twist-4 four guark operator to be, 

O^u = (gr^g)(gr^g) - ^gf,u{q^aq)iq^aq), (21) 

then for (p = 0|(9oo|p = 0), the matrix element can be written as the sum of full sguares: 

{p = 0\[l{qToqy + lmqnp = 0) (22) 

Hence both M4 and A should be positive. Since both come from four guark operators, we 
expect them to be similar in magnitude. We will vary M4 from to 2 x A and solve for A 
and B. The result is summarized in Table 1. 
Our analysis suggest that: 

1. The magnitude of the guark gluon mixed operator B is relatively large compared to 
the four guark operators at the 5 GeV^ scale. The former is in the range of -(350~450 
MeV)^ and the latter (100 ~ 200 MeV)^, consistent with previous estimates 0. 

2. C is negligible compared to B but not so compared to A. For the transverse moment 
in eg.(|l9D, the Wilson coefficient of B is small and its contribution is largely cancelled 
by that of A. Therefore, the contribution from C, which comes with a large Wilson 
coefficient, is a sizable part of the estimated twist-4 contribution and can not be 
neglected. 



IV. CONCLUSION 

We have considered the contribution of guark-mass-dependent twist-4 operators to the 
Bjorken sum rule and momentum sum rule. In the Bjorken sum rule only operators propor- 
tional to appear. This is easily seen because the operator linear in m would be qrri'j^Daq 
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which vanishes when using the equations of motion. So the correction to the leading part is 

2 

^ and negUgible. 

However, for the momentum sum rules, qDaDpmq contributes with a large Wilson coef- 
ficient and becomes a non-negligible correction to other twist-4 effects. 
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Appendix A 

In this appendix, we summarize the form of quark propagator in external gauge field 
including mass dependent terms relevant for the moment calculations. We will group them 
together according to the dimension of external gauge field plus the dimension of quark 



mass. 



Such that 



S{p) = + + 5(p)(2) + ^(p)(3)_ (A.l) 



p2 



= Lp^ 7/375 + \m^Pa7a, 

2 1 

6 p 

-Sip^D^Pa Gap 7/375] - -^rngGapcTap (A. 2) 

Zp 



Appendix B 



In this appendix, we summarize useful identities that are needed to reduce eq.(|TO[) into 
its finial form with independent matrix elements. We will only look at the spin-1 part of 
the following operators. 

q DaDfsD^ -fsq = Qaf^W^D - A - 3B - 2C].,s + 9ajU^D + 5A - 3B - 2C]ps 
+ QasU-'^D -A + B + QC]p^ + gp^U2D -A-W- 2CU 
+ gp&W^D + A + B + QCU + g^^U^'^D -A + B + QC]pa, (B.l) 

where we have neglected the terms proportional to e tensor. 
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q DaD/3D^ isl^q = e^^^5^[-25 - A - D + AC]^^ + e«,^5^[-2A - 2D + 8C]^/3 

+ e^p^sU-'^B -A-D + 4C]^., (B.2) 



where we have neglected the terms proportional to g. Using eg. ( |B.2|) , we can show, 

E = q{iD„ ^.4757/3? =[-B-A-D + AC]^p (B.3) 

Other straightforward expansion are, 

< — ^ -1 -1 

+ gp&^^[W + E]ps - g^&U^B + EU, (B.4) 

and, 

qDaGisy'ysq = fi'a/3|[-3A - D]^s - 5(„^|[-3yl - D]ps 

+ g(ssl[3D + A]^s - g^sH^D + A^. (B.5) 
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M4 


A 


B 





0.041 GeV2 


-0.148 GeV2 


A = 0.017 GeV^ 


0.017 GeV^ 


-0.165 GeV^ 


2A = 0.022 GeV^ 


0.011 GeV^ 


-0.169 GeV2 



Table 1: Proton matrix elements of M4, A and B. 
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